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Abstract. We present a power spectrum analysis of the long ASCA observation of Ark 564 in June/July 2001. 
The observed power spectrum covers a frequency range of ~ 3.5 decades. We detect a high frequency break at 
~ 2 X 10“® Hz. The power spectrum has an rms of ~ 30% and a slope of ~ —1 and ~ —2 below and above 
the break frequency. When combined with the results from a long RXTE observation (Pounds et al. 2001), the 
observed power spectra of Ark 564 and Cyg X-1 (in the low/hard state) are almost identical, showing a similar 
shape and rms amplitude. However, the ratio of the high frequency breaks is very small (~ 10^“'*), implying that 
these characteristic frequencies are not indicative of the black hole mass. This result supports the idea of a small 
black hole mass/high accretion rate in Ark 564 . 
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1. Introduction 

Narrow-line Seyfert 1 (NLSl) galaxies are a peculiar group 
of active galactic nuclei (AGN) characterized by their dis¬ 
tinct optical line properties (Osterbrock & Pogge 1985). 
In hard X-ray studies NLSl galaxies comprise less than 
10% of the Seyfert galaxies, however, from the ROSAT 
All-Sky Survey it became clear that about half of the 
AGN in soft X-ray selected samples are NLSl galaxies 
(Grupe 1996, Hasinger 1997). Boiler et al. (1996) found 
from ROSAT observations that the soft X-ray spectra of 
NLSl galaxies are systematically steeper than those of 
broad line Seyfert 1 galaxies. They further discovered that 
NLSl galaxies frequently show rapid short time scale X- 
ray variability which can be interpreted as evidence for a 
small black hole masses in these objects. 

Ark 564 is the X-ray brightest NLSl galaxy with a 
2—10 keV flux of ~ 2 X 10~^^ ergcm“^ s“^(Vaughan et al. 
1999) and shows large amplitude variations on short time 
scales (Leighly 1999). Therefore it is the best candidate 
to study its X-ray variability in order to obtain impor¬ 
tant clues about the size of the black hole mass and the 
accretion rate in NLSl galaxies. Recently, Ark 564 was 
observed by RXTE once every ~ 4 days from January 
1999 to September 2000 covering ~ 20 months of data. Its 


2 — 10 keV power spectral density (PSD) function showed 
a cut-off at a frequency which corresponds to a time scale 
of ~ 13 days (Pounds et al. 2001). Ark 564 was also ob¬ 
served for a period of ~ 35 days in June/July 2000 by 
ASCA . This observation was part of a multi-wavelength 
AGN Watch monitoring campaign (Turner et al. 2001). 
The flux variations that Ark 564 exhibited during this 
observation have already been studied by Edelson et al. 
(2001) who found that the variability amplitude is almost 
independent of energy band and the power spectrum is 
harder in the hardest energy bands. There were no de¬ 
lays between the variations in the different bands and no 
signs of non-linear behaviour in the light curves. Gliozzi et 
al. (2001) also found no sign of non-linear behaviour and 
no statistically significant indication of non-stationarity in 
the light curves. Furthermore, using nonlinear techniques 
they were able to demonstrate that the source behaves 
differently in the high and low flux states. 

In this work we present a timing analysis of the 
June-July ASCA observation of Ark 564 based on power 
spectral analysis techniques. In the following section we 
present the data and the results from the PSD analysis 
while in Section 3 we discuss our conclusions. 
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Fig. 1. Plot of the ASCA SISO+1 and GIS2+3 0.7-, 
keV light curves (open and filled circles, respectively). Th 
data are binned in 5400 sec intervals (~ the ASCA orbita 
period). For clarity, the SIS light curve is shifted by +: 
counts/sec along the y—axis. 

2. Data Analysis 

ASCA observed Ark 564 from 2000 June 01 11:51:27 ti 
2000 July 05 23:57:54 (UT). We extracted the X-ray data 
from the public archive at ISAS and applied standard cri¬ 
teria for the data analysis in order to create light curves for 
all four instruments. Details of the data reduction process 
are given in Brinkmann et al. (2001). 

Light curves were extracted in the energy bands 0.7- 
5 keV, 0.7-2 keV (soft band), 2-5 keV (medium band) 
and 5-10 keV (hard band) for both instruments. In all 
cases we combined the data from the SIS and CIS detector 
pairs. Figure 1 shows the SISO+1 and CIS2+3 0.7-5 keV 
light curves in 5400 s bins (~ the ASCA orbital period). 
Clearly, the source is significantly variable. The root mean 
square variability amplitude (corrected for the experimen¬ 
tal noise) is 33%, 32% and 29% for the soft, medium and 
hard band light curves, respectively. 

2.1. Power Spectrum Analysis 

We used the 0.7-5 keV CIS and SIS light curves (which 
have the highest signal to noise ratio) to estimate the over¬ 
all X-ray PSD of Ark 564 . First, we used the 16 sec binned 
SIS and CIS light curves to compute the high frequency 
part of the PSD (i.e. the PSD at frequencies > 10“^’® Hz). 
There are 476 and 531 segments in the SIS and CIS light 
curves, respectively, with no gaps in them. Their duration 
is between ~ 1000 and ~ 3500 sec. For each part we com¬ 
puted the periodogram and normalised it to the square 
of the mean count rate. Then, we combined all the peri- 
odograms in one file, we sorted them in order of increasing 
frequency, we computed their logarithm and grouped them 
into bins of size 500 following the method of Papadakis & 



Fig. 2. Plot of the 0.7-5 keV SISO+1 and CIS2+3 power 
spectra (open and filled circles respectively). The dotted 
and dashed lines show the best fitting power law model 
lines to the SIS and CIS PSDs, respectively. The solid 
lines show the expected Poisson noise power level for the 
high frequency PSDs. 

Lawrence (1993). To compute the PSD at low frequencies 
(i.e. < 10“^ Hz) we used the 5400 sec binned light curves. 
The use of the approximate ASCA orbital period as the 
bin size results in an evenly sampled light curve with very 
few missing points (10 and 8 out of 552 in the CIS and 
SIS light curves, respectively). We accounted for them us¬ 
ing linear interpolation between the bins adjacent to the 
gaps. We computed the periodograms, and averaged their 
logarithm using a bin size of 20. 

The low and high frequency PSDs of the SIS and CIS 
light curves are shown in Figure 2. The SIS and CIS power 
spectra look similar. Using standard statistics we fitted 
a power law model of the form log[P(:z)] = log[Aj^““+C'„], 
separately to the high and low frequency part of the PSDs 
(we denote with a\i and Ohf the slope of the low and high 
frequency parts respectively). The constant Cn represents 
the contribution of the experimental Poisson noise process 
to the observed PSD. It can be estimated using the aver¬ 
age error of the points in each light curve (e.g. Papadakis 
& Lawrence 1995) and, as a result, Cn was kept fixed at 
its expected value (shown with solid lines in Figure 2 for 
the high frequency power spectra) during the model fit¬ 
ting procedure. The best fitting results are listed in Table 
1 and the best fitting models are plotted in Figure 2 (dot¬ 
ted and dashed lines for the SIS and GIS power spectra, 
respectively). All errors quoted in this paper correspond 
to the 68% confidence region and were computed following 
the method of Lampton et al. (1976). 

The best fitting slope values of the SIS and GIS power 
spectra are consistent within their errors (a®P — ap/® = 
-0.15 ± 0.07 and af/® - aff^® = 0.08 ± 0.18). Therefore, 
despite the slight difference in the responses of the detec¬ 
tors on board ASCA , the power spectra of the GIS and 
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Table 1. Results of the power law model fits to the 0.7-1 
keV, SIS and GIS power spectra of Ark 564 . 



ait 

Ohf 

Xff/dof 

Xhf/dof 

SISO+1 

1.31 ±0.13 

1.65 ±0.04 

15.6/11 

78/57 

GIS2-t3 

1.23 ±0.13 

1.80 ±0.06 

13.2/11 

106/59 


SIS light curves are indeed similar. The weighted meai 
Ohf and a\i values are 1.70 ±0.03 and 1.27 ±0.09, respec 
lively. Their difference, 0.43 ± 0.10, shows that the lo'v 
and high frequency parts of the power spectra have sig 
nificantly different slopes. As Figure 2 shows, the slope o 
the PSD changes at ~ 10“^ Hz. Above that frequency 
the slope is significantly steeper than the slope at lowe 
frequencies. 

Since the SIS and GIS power spectra are similar, w 
combined them to produce the final, 0.7-5 keV PSD o 
Ark 564 . The resulting power spectrum (with the Poissoi 
noise subtracted) is shown in Figure 3. We have groupec 
the combined GIS and SIS periodograms into bins of sizi 
200, 600 and 1000 for the high frequency part of the spec 
trum (the bin size is increased towards the higher frequen 
cies in order to increase the signal to noise) and into bin 
of size 50 for the low frequency part. A simple power la^; 
model cannot fit the PSD well (y^ = 259 for 58 dof) 
For that reason, we used a “broken power law” model 
According to this model, the low frequency part of the 
PSD is fit to a function of the form and, above 

a “break frequency” the PSD is fitted to a power law 
with a different index, Ohigh- The fit results are listed in 
Table 2; the best fit model is plotted in Figure 3 (solid 
line). Gompared to the simple power law model, a broken 
power law provides a significantly better fit to the data. 
The PSD of Ark 564 has a ~ —1.2 slope at frequencies up 
to ~ 2 X 10“^ Hz above which the PSD steepens to a slope 
of ~ —2. 

Finally, we investigated whether the PSD parameters 
depend on energy. We used the 0.7-2 keV, 2-5 keV and 
5-10 keV light curves to compute the power spectra in 
the three bands. A power law model cannot fit the soft 
and medium band PSDs (y^ = 206 and 157 for 66 and 
53 dof, respectively). The broken power law model gives 
a significantly better fit to these power spectra. The hard 
band PSD cannot be computed for frequencies higher than 
a few X10“^ Hz due to the low signal to noise of the 16- 
sec hard band light curves; therefore the simple power law 
model gives already an acceptable fit to the PSD in this 
case. The best fitting results for the power law (for the 
hard band PSD) and broken power law models (for the 
soft and medium band PSDs) are listed in Table 2. The 
model parameters of the soft and medium energy PSDs are 
consistent within the errors, except from oiow Our results 
show that the low frequency PSDs become significantly 
harder with increasing energy. This is in agreement with 
the results of Edelson et al. (2001). 
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Fig. 3. Plot of the combined GIS and SIS 0.7-5 keV power 
spectra. The solid line shows the best fitting “broken 
power law” model to the data. The dashed line is the 
best fitting model to the RXTE PSD (from Pounds et al. 
2001 ). 


Table 2. Results of the “broken power law” model fits to 
the PSDs of Ark 564 


Energy 

Qllow 

Q^high 

Vht 

X^/dof 




(xl0“® Hz) 


0.7-5 keV 
0.7-2 keV 
2-5 keV 
5-10 keV 

1 

-‘-•^^- 0.04 

-l-.^0_o.o4 

1 1 q+0.03 

J---Lo_o.o 4 
-1 -, 0 + 0.06 
J---LO o.o 3 

“^“qFTTTTT 

l.yU_o.i 4 

-L-c? 0 _q 20 

1 07 + 0.34 
-L-0/_o.32 

Q q + U.b 

^•o_ 0.6 

1 '7+0.6 
^•'- 0.9 

72.6/56 

84.1/64 

84.1/51 

52.9/31 


3. Discussion 

Using the ASCA June-July 2001 observation of Ark 564 
we have estimated its power spectrum over a frequency 
range of ~ 3.5 decades. We have detected a high frequency 
break at ^ 2 x 10“^ Hz. 

In Figure 3, together with the 0.7—5 keV PSD, we have 
also plotted the best fitting “cut-off power law” model 
from Pounds et al. (2001). The observed PSD of Ark 564 
extends over 6 decades in frequency. It is clear from 
this Figure that the PSD starts to flatten below ^ 10“® 
Hz, showing a second, low frequency break. MGG 6-30- 
15 also shows two break frequencies in its PSD (Nowak & 
Ghiang 2000), however, the evidence for the low frequency 
break is not significant in that case (Uttley et al. 2001). 
Therefore, Ark 564 is probably the only Seyfert galaxy 
for which we are certain that the observed PSD has a 
shape that is almost identical to the shape of the Gyg X-1 
PSD in the low/hard state. In both cases, the PSD slope 
is fiat below a characteristic low frequency. Above this 
frequency the PSD steepens to oc and then becomes oc 
v~^ above a second break frequency. Apart from the PSD 
shape, the rms amplitudes at frequencies higher than the 
low frequency cut off are also very similar. These results 
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show that Ark 564 corresponds to the low/hard (instead of 
the high/soft) state of the Galactic black hole candidates 
(GBHC). 

However, despite this similarity, Ark 564 is not simply 
a “larger” version of GBHCs. The two break frequencies in 
Cyg X-1 are located between 0.03 — 0.3 and 1 — 10 Hz (e.g. 
Nowak et al. 1999, Negoro et al. 2001). The ratio between 
the characteristic frequencies in Ark 564 and Cyg X-1 is 
therefore ~ 10®“® and ~ 10®“^, if we consider the low 
and high frequency breaks, respectively. Although there is 
an uncertainty associated with these ratios (for example, 
we do not know whether the break frequencies in AGN 
are constant or change with time) this cannot probably 
explain the two orders of magnitude difference between 
the ratios. Perhaps then the break frequencies do not scale 
proportionally with the black hole mass. For example, if 
the high frequency break in Cyg X-1 corresponds to the 
Keplerian orbital period at ~ 30i?s around the central 
object (Nowak et al. 1999; Rs is the Schwarzschild radius) 
and the high frequency break in Ark 564 corresponds to 
the same time scale but at a radius ~ 3i?s, then the ratio 
of the characteristic time scales will not be equal to the 
mass ratio but to the mass ratio x (30i?s/3i?s)“®/^. 

On the other hand, if the low frequency breaks scale 
with mass, then the central black hole mass of Ark 564 
should be ~ 10^ Mq, assuming that the black hole in Cyg 
X-1 is 10 Mq (Herrero et al. 1995). If that is the case, 
taking into account the luminosity of Ark 564 , Pounds 
et al. (2001) concluded that the source is accreting at a 
substantial fraction of the Eddington limit. The orbital 
period at 3Rs for a 10^ Mq black hole is ~ 5000 sec, 
which is an order of magnitude larger than the time scale 
which corresponds to the high frequency break in the PSD 
of Ark 564 . Either the mass of the black hole in Ark 564 
is smaller than 10^ Mq (in which case the accretion rate 
is larger than the Eddington limit), or due to the large 
accretion rate, the characteristic time scales are smaller. 
When the accretion rate approaches the Eddington limit 
advective energy transport dominates over radiative cool¬ 
ing and the disk becomes moderately geometrically thick, 
a so called “slim disk” (Abramowicz et al. 1988). In this 
case, a substantial amount of radiation can be produced 
from inside the last stable orbit, i.e. 3Rs (Mineshige et 
al. 2000). The free-fall time scale at distance R is given 
by tfi = {Rg/c){R/Rg)^^^ = 500Ms{R/R g)^^^ sec (Rees 
1984), where Mg is the mass in units of 10®Mq and 
Rg = GM/c^ is the gravitational radius. For Mg = 0.1 
and R = 3Rs = 6 Rg, we get tg ~ 700 sec which is compa¬ 
rable to 1/ i%f that we find for Ark 564 . The sound cross¬ 
ing time scale is another possible candidate that could cor¬ 
respond to l/ri 3 f. In the thin disk case this time scale is 
very long even in the innermost parts of the disk. However, 
when the accretion rate is comparable to the Eddington 
rate and the disk is dominated by radiation pressure, the 
group velocity of sound waves (responsible for transmit¬ 
ting density fluctuation information) becomes very large 
(e.g. Krolik et al. 1991). Gonsequently, the sound crossing 


time scale at say R = 3Rs will be decreased substantially 
and could be of the order of l/i^f. 

We conclude that, the similarity of the PSD and time 
lags in Ark 564 implies that X-rays in this source are pro¬ 
duced by a mechanism similar to the X-ray emission mech¬ 
anism in Cyg X-1 in its low/hard state. On the other hand, 
the comparison between the characteristic time scales in 
the two systems suggests that at least one physical pa¬ 
rameter is different in them, probably the accretion rate, 
with Ark 564 having a substantially higher accretion rate 
than Cyg X-1 in its low/hard state. 
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